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ABSTRACT 

Galaxies and quasars are thought to provide the bulk of the photons responsible 
for ionizing the hydrogen in the intergalactic medium (IGM). We use a large set of 
hydrodynamical simulations, combined with measurements of the Lya opacity of the 
IGM taken from the literature, to obtain robust estimates of the photoionization rate 
per hydrogen atom at z — 5 and 6. We find the photoionization rate drops by a factor 
of two and four, respectively, compared to our recent measurements at z — 2 — 4. The 
number of ionizing photons emitted by known sources at z — 5 and 6, based on an 
extrapolation of source numbers below the detection limit and standard assumptions 
for the relationship between the ionizing emissivity and observed luminosity density 
at 1500A, are in reasonable agreement with the photoionization rates inferred from 
the Lya forest if the escape fraction of ionizing photons from galaxies is large ( > 20 
per cent). The expected number of ionizing photons from observed sources at these 
redshifts therefore appears sufficient to maintain the IGM in its highly ionized state. 
Claims to the contrary may be attributed to the adoption of an unduly high value for 
the clumping factor of ionized hydrogen. Using physically motivated assumptions for 
the mean free path of ionizing photons our measurements of the photoionization rate 
can be turned into an estimate of the ionizing emissivity. In comoving units the inferred 
ionizing emissivity is nearly constant over the redshift range 2 — 6 and corresponds 
to 1.5 — 3 photons emitted per hydrogen atom over a time interval corresponding to 
the age of the Universe at z = 6. This strongly suggests that the epoch of reionization 
was photon-starved and extended. Completion of reionization at or before z = 6 
requires either an emissivity which rises towards higher redshifts or one which remains 
constant but is dominated by sources with a rather hard spectral index. For standard 
assumptions, the ionizing emissivity required for completion of reionization at or before 
z — 6 lies at the upper end of recently reported values from searches for high redshift 
galaxies at z — 8 — 10. 

Key words: hydrodynamics - methods: numerical - intergalactic medium - quasars: 
absorption lines - diffuse radiation. 



1 INTRODUCTION 

Hydrogen in the intergalactic medium (IGM) is maintained 
in its highly ionized post-reionization state by the meta- 
galactic ionizing background. This radiation is attributed to 
the integrated UV emission shortward of the Lyman limit 
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from young star forming galaxies and quasars, which is sub- 
sequently filtered and reprocessed as it propagates through 
the clumpy IGM (lHaardt fc Madaulll99rJ : iFardal et all 19981 : 
iMadau et afl ll999h Observational constraints on the am- 
plitude and spectral shape of the ionizing background 
can be used to infer the relative abundances of differ- 
ent UV photon emitting populat ions and their evolution 
(iMiralda-Escude fc Ostrikerlll990l; iMeiksin fc Madad ll99~' 



iBi fc Davidsenl Il997l ; iDevriendt et alj Il998l : IShull etal 
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1999: Bianchi et all 2001; Haehnelt etafl l200ll; IShull et all 



2004 



Demiariski fc Doroshkevichl |2004| ; iBolton et alj 120051 . 
20061 ). The integrated UV emissivity from galaxies can also 
be directly linked to the star formation density i n the 
Universe at a particular epoch (jMadau et al.l Il998). and 
the metagalactic hydrogen ionization rate is closely related 
to the evolution of the neutral hydrogen fraction in the 
IGM (|Fan et al.ll2002l ; iFan et al. 1120061 ; feecker et alj|2007l ). 
Models of the metagalactic ionizing background are also 
an essential ingredient w ithin hydrodynamical simulations 
of st r ucture format i on JCenl 19921; Miralda-Escude et alj 
1 19961 : iTheuns et all Il998l : I Jena et al] l2005h and provide 
the ionization corrections needed to infer the abundances 
of heavy ele ments in the IGM f r om metal ion ab- 
sorption lines llGiroux fc Shulll 1 19971 : lAguirre et aD 120051 ; 
lOppenheimer fc Davell2006h . 

Consequently, much effort has been directed towards 
constraining the nature and evolution of the ionizing back- 
ground over a wide range of redshifts. Numerical simula- 
tions of structure formation, combined with the Lya opac- 
ity distribution observed in quasar absorption spectra, 
have enabled significant progress to be made in this 
field dRauch et al.fl 19971; iMcDonald fc Miralda-Escudell200ll; 



of the Lya opacity in high redshift quasar spectra (|Songailal 
l2004l : lFan et al. 1120061 ) , combined with a large suite of high- 



Meiksm fc White! 12003: iTvtler et all 120041: IBolton et aD 



20051 : Jena et all 120051 ). However, at redshifts approaching 



and beyond the putative tail-end of the hydrogen reioniza- 
tion epoch at z ~ 6, determining the ionization s tate of 
the IGM is still problematic ( e. a. IFan et all 120021: ISantosi 
20041; ISongailall2004|: iGallerani et a.l.ll2006l: IFan et al. 1120061 : 



Bolton fc HaehneltlboOTl ; iDiikstra et al.l 120071 ). It is gener- 
ally agreed that by z ~ 6 the ionizing radiation emit- 
ted by quasars alone is insufficient to reionize the IGM 
llFan et al]|200ld;lMiralda-Escudell2003l;ISchirber fc Bullockl 
2003 1; IDiikstra et al]|2004l; lYan fc Windhorsdl2004 iMeiksinl 
2005l ; ISrbinovskv fc Wvithell2007l ; IShankar fc Mathurll2007l ). 
In addition, there have been some claims that the ob- 
served galaxy population at 2 ~ 6 is also unable 
to produce the number of hydrogen ionizing photons 
needed to maintain the surrounding IGM in a highly ion - 
ized state (|Lehnert fc Bremer! 120031 ; iBunker et all l2004h . 
However, other studies have i ndicated this deficit may 
not real dStiavelli et aL I |2004| ; lYan fc Windhorstl [2004b; 
iBouwens et all 120061 ) . and may instead be attributable to 
an overestimate of the number of ionizing photons required 
to balance recombinations in the IGM. There may also be 
a population of fainter galaxies which remain undetected. 
Considerable uncertainties are also introduced when con- 
verting the galaxy luminosity function into a UV emissivity 
and hence an ionizing photon production rate. An assump- 
tion must also be made for the very uncertain escape frac- 
tion of io nizing photons from their respective galactic envi- 
ronments llGiallongo et al]|2002l; iFernandez-Soto et al ]|2003l ; 
llnoue et al. 1 120061; ISiana et al. 1120071 ') , which is likely to vary 
considerably (|Shaplev et al]|2006h . 

In this paper we shall readdress the question of whether 
or not the observed ionization state of the IGM at z ~ 
5 — 6 is consistent with the ionizing emissivity inferred 
from the observed space density of star-forming galaxies 
and quasars. We approach this issue with improved mea- 
surements of the metagalactic photoionization rate per hy- 
drogen atom, T_i2 = Thi/10 -12 s _1 , at z — 5 and 6. 
These are obtained from recently published measurements 



resolution hydrodynamical simulations (e.g. iRauch et al 
19971; IMcDonald fc Miralda-Escude 200ll;ICen fc McDonald 



20021 ; ITvtler et al] |2004| ; IBolton et al] 120051 )7 These esti- 
mates represent the most detailed constraints on T_i2 ob- 
tained at these redshifts using state-of-the-art hydrodynam- 
ical simulations, and the systematic uncertainties are care- 
fully explored for a wide range of model parameters. These 
new data are then compared to the ionization rate inferred 
from our own estimates for the ionizing emissivity derived 
from published measurements of the space density of Lyman 
break galaxies (LBGs) and quasars, combined with a sim- 
ple model for the mean free path of ionizing photons. We 
also c ompare this approach to that used by iMadau et al.l 
(1999), who instead compared an estimate of the observed 
ionizing emissivity to that required to balance recombi- 
nations in a simple model where the inhomogeneous spa- 
tial distribution of the IGM is characterised by an overall 
H 11 clumping factor. In this way, we derive an upper limit 
on the H 11 clumping factor from the T_i2 we measure at 
2 = 6. Finally, we turn our new measurements of the ion- 
ization rate at z ~ 5 — 6 into a measurement of the ionizing 
emissivity. Given this constraint, we then explore the range 
of plausible evolutionary histories for the ionization state 
of the IGM towards yet higher redshifts and compare these 
to recent observational estimates of the UV emissivity at 
z = 8 - 10. 

The structure of this paper is as follows. In Section 2 
we describe our procedure for estimating the metagalactic 
hydrogen ionization rate at z — 5 and 6 using measurements 
of the IGM Lya opacity and hydrodynamical simulations of 
the IGM. We pay particular attention to the systematic un- 
certainties involved in such a measurement, an d as such this 
work closely follows the methodology used by IBolton et al] 
(2005) (hereafter B05). We refer the reader to this work for 
further details. In Section 3 we discuss the model we use 
for the mean free path of ionizing photons, and in Section 
4 we use published LBG and quasar luminosity functions to 
directly compute the Lyman limit emissivity at z = 5 and 
6. We discuss a model for spatial fluctuations in the ioniz- 
ing background in Section 5. This is used to estimate the 
possible systematic effect these have on our measurements 
of r_i2. In Section 6 we compare the metagalactic hydro- 
gen ionization rates inferred from both the Lya opacity of 
the IGM and the combined galaxy and quasar Lyman limit 
emissivities, and in Section 7 we discuss the correspond- 
ing comoving ionizing emissivity and its implications for the 
ionization state of the IGM at z > 6. Implications for fu- 
ture 21cm experiments and the evaporation of mini-haloes 
are also considered. Lastly, in Section 8, we summarise and 
present our conclusions. 



2 INFERRING THE IONIZATION RATE 
FROM THE IGM Lya OPACITY 

2.1 Simulations of the IGM 

The hydrodynamical simulations used for this st udy were 
run using the parallel TreeSPH code GADGET- 2 (|Springell 
2005). Our fiducial simulation volume is a periodic box 
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Table 1. Simulations used in our study of the dependence of 
r_i2 on various cosmological and astrophysical parameters. A 
flat Universe with = 1 — f2 m is assumed. The last column lists 
the factor by which the He II photoheating rate is multiplied to 
investigate the effect of gas temperature on the inferred T_i2. All 
simulations listed have a box size of 15ft - 1 comoving Mpc and 
contain 2 X 200 3 gas and dark matter particles. 



Name 




tt b h 2 


ft 


<T8 


n 


AhcII 


15-200 


0.26 


0.024 


0.72 


0.85 


0.95 


1 


Tl 


0.26 


0.024 


0.72 


0.85 


0.95 


0.2 


T2 


0.26 


0.024 


0.72 


0.85 


0.95 


0.5 


T3 


0.26 


0.024 


0.72 


0.85 


0.95 


1.5 


T4 


0.26 


0.024 


0.72 


0.85 


0.95 


2 


Ml 


0.17 


0.024 


0.72 


0.85 


0.95 


1 


M2 


0.40 


0.024 


0.72 


0.85 


0.95 


1 


M3 


0.70 


0.024 


0.72 


0.85 


0.95 


1 


M4 


1.00 


0.024 


0.72 


0.85 


0.95 


1 


SI 


0.26 


0.024 


0.72 


0.50 


0.95 


1 


S2 


0.26 


0.024 


0.72 


0.70 


0.95 


1 


S3 


0.26 


0.024 


0.72 


1.00 


0.95 


1 


Si 


0.26 


0.024 


0.72 


1.20 


0.95 


1 



15/i -1 comoving Mpc in length containing 2 x 200 3 gas 
and dark matter particles. We find this provides the best 
compromise between accuracy and speed. Star formation is 
included using a simplified prescription which converts all 
gas particles with overdensity A > 10 and temperature 
T < 10 s K into collisionless stars . The uniform UV back- 
ground model of lHaardt fc Madaul l|200ll ) (hereafter HM01) 
including emission from quasars and galaxies is employed in 
the optically thin limit, and the ionization state of the gas 
is computed using the non-equilibrium ionization algorithm 
of lBolton et al.l l|2006t ). The simulations were all started at 
z — 99, wit h initial conditions gene rated using the transfer 
function of lEisenstein fc Hul (|l999h . The cosmological pa- 
rameters of the simulations used in this study are listed in 
Table 1, with the final column listing the factor by which the 
He ii photo-heating rate is artificially raised to explore the 
effect of gas temperature on the inferred F_i2. The cosmo- 
logical parameters used for the fiducial model are consistent 
with the combin ed analysis of the third year WM AP and 
Lya forest data (|Viel et al.ll2006l ; ISeliak et alj|2006l '). 

To check numerical convergence we also run six fur- 
ther hydrodynamical simulations with differing box sizes 
and mass resolutions. The different resolution parameters 
for these simulations are listed in Table 2. All other aspects 
of these simulations are identical to the fiducial 15 — 200 run. 



2.2 Synthetic spectra generation and fiducial 
parameter ranges 

Synthetic Lya spectra are constructed at z — 5 and 6 us- 
ing 1024 random lines-of-sight drawn parall el to the box 
boun daries from each of the simulations (e.g. iTheuns et all 
1998). Each line-of-sight consists of 1024 pixels. To deter- 
mine F_i2 we rescale the synthetic spectra to match the 
mean normalised flux, (F), observed in the Lya forest por- 
tion of quasar spectra at z = 5 and 6. This is achieved by 
linearly rescaling the optical depths in each pixel of the syn- 
thetic spectra by a constant factor, A, such that 



Table 2. Resolution and box size of our six additional simulations 
which have the same parameters as the 15-200 model used for 
our parameter study. The mass resolution used for our parameter 
study simulations is also listed for comparison. 



Name 


Box size 


Total particle 


Gas particle 




comoving Mpc 


number 


mass [/i _1 Mq] 


15-400 


15ft" 1 


2 X 400 3 


6.78 x 10 5 


30-400 


30k- 1 


2 X 400 3 


5.42 x 10 e 


15-100 


15ft- 1 


2 X 100 3 


4.34 x 10 7 


30-200 


30ft" 1 


2 X 200 3 


4.34 x 10 7 


60-400 


60ft- 1 


2 x 400 3 


4.34 x 10 7 


30-100 


30ft- 1 


2 x 100 3 


3.47 x 10 s 


15-200 


15ft" 1 


2 x 200 3 


5.42 x 10 6 



10.0 



OFon et ol. (2006) 
ASongoilo (2004) 
□ Kirkmon et ol. (2005) 
• Schoye et ol. (2003) 
▼ This work 



1.0 




0.1 DD 



Figure 1. Observational constraints on the Lya effective opti- 
cal depth of the IGM and its evolution with redshift. The open 
squares, triangles and diamonds correspond to the data of Kirk- 
man et al. (2005), Songaila (2004) and Fan et al. (2006) respec- 
tively. The filled circles with error bars show the estimates for T ef r 
used in B05, based on the data of Schaye et al. (2003). The filled 
inverted triangles with error bars correspond to the constraints 
on T e ff we use for this work. 



JV 



(1) 



where 73 is the optical depth in each of the N pixels in the 
synthetic spectra and r e ff = — ln(_F) is the observed Lya ef- 
fective optical depth. The F_i2 required to reproduce r e g is 
then given by T_i2 = T a i m /A, where r s i m is the ionization 
rate originally used in the simulation. Repeating this pro- 
cedure for many simulations with varying cosmological and 
astrophysical parameters allows one to d etermine how r_i 2 
scales with these parameters (e.g. B05. [jena et al.ll2005l ). 
The values adopted for r e s in this st u dy ar e base d on the 
recent data published by iFan et al~l ()2006l ) and ISongailal 
(2004). These data a re displayed in Fig. [T] along with the 
lower redshift data of iKirkman et all ((20051). The filled cir- 
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Table 3. Fiducial parameter values and estimates of their uncer- 
tainties. The values listed for r e g are at z = [5,6]. 

Parameter Fiducial values and uncertainties 



Table 4. The redshift dependent indices for the scaling relations 
between T_i2 and several cosmological and astrophysical param- 
eters. 



T 

T e S 

7 

tt h h 2 
is 
h 



1 ±0.5 x 10 4 K 
0.26±0.04 
[2.07to'j?.> 5-50] 
1.3±0.3 
0.024±0.001 
0.85±0.05 
0.72±0.04 



cles correspond to T e s measured by ISchave et al.l (|2003h . 
used in B05 for determining T_i2 at 2 < z < 4. 

Approaching z = 6 there is some debate whet her 
an abrupt transition (IFan et al.l 120 02; Fan et al. 20061) or 
smooth progression (ISongailal 12004 ; iBecker et al.l 120071 ) in 
the redshift evolution of r c ff is a better fit to the observa- 
tional data. Rather than attempt to parameterise the red- 
shift evolution of r c ff, we bin the data into redshift bins 
of width Az = 0.25 and compute the mean in each bin. 
These are shown at z — 5 and 6 as the filled inverted tri- 
angles in Fig. [T] with corresponding uncertainties estimated 
from the interquartile range of the binned data. Note that at 
z ~ 6 some quasar sight-lines exhibit a full lGunn fc Peterson! 
(1965) trough, and therefore only a lower limit on T eff , and 
hence upper limit on F_i2, may be obtained (|Fan et al. I 
2006). The fiducial values for the Lye* effective optical depth 
adopted for this study are therefore T e g = 2.07^Q27 at z — 5 
and a lower limit of T e « > 5.50 at z — 6, corresponding to 
(F) = 0.127lo;°27 and (F) < 0.004, respectively. 

In addition, we must also make some assumptions 
for the other simulation input parameters on which the 
IGM Lya opacity depends. The adopted fiducial values and 
uncertainties for the cosmological parameters are S7 m = 
0.26 ± 0.04, Q h h 2 = 0.024 ± 0.01, h = 0.72 ± 0.04 and 
as = 0.85 ± 0.05. These are consistent with the combined 
analysis of the t hird year WMAP and Lya forest data 
l|Viel et all |2006| ; ISeliak et al.l 120061 ). Unfortunately, there 
are currently no reliable constraints on the thermal state 
of the IGM at 2 > 4. However, the slope of the effective 
equation of state for the low density IGM, T — To A 7-1 , 
is likely to be well bracketed by assumin g 7 = 1.3 ± 0.3 
l|Hui fc Gnedinlll997l ; IValageas et"aHl2002ft . We adopt this 
as the fiducial range for 7 in this study. The temperature 
of the IGM at mean density, To, is more difficult to place 
limits upon. We shall assume that the double reionization of 
helium does not occur at the same time as hydrogen reion- 
ization; there is evidence for the tail end of He 11 reionization 
occurring around z p± 3 (ISongaila Il99£ ; ISchave et al. 2000; 
iBernardi et~aLll2003l ; iBolton et al.ll2006l ). Therefore, T may 
reach 1.5 — 2 x 10 4 K if hydrogen a nd single helium reion - 
ization occurred just above 2 = 6 (|Hui fc Haimanl [2003). 
although if hydrogen reionization occurred earlier To may 
be as low as 0.5 x 10 4 K. We adopt a fiducial range of 
1 ± 0.5 x 10 4 K for this study. All our fiducial parameters 
and their uncertainties are summarised in Table 3. 



z 


To 


7 


S7 m 


0"8 Toff 


5 


-0.57 


0.81 


-1.27 


-1.63 -1.74 


6 


-0.61 


1.07 


-1.38 


-1.90 -1.88 



2.3 Numerical resolution 

We must also assess the extent to which changes in the nu- 
merical parameters of the simulations, in this case box size 
and resolution, affect the F_i2 inferred in our study. This 
is achieved by using the simulations with different box sizes 
and mass resolutions listed in Table 2. 

Comparing T_i2 inferred from the 15 — 100 and 60 — 400 
simulations, which have the same mass resolution, the larger 
box size lowers T_i2 by 10 per cent at 2 = 5 and 12 per cent 
at 2 = 6. A smaller reduction of 4 per cent at both redshifts 
is found on comparing the 30 — 200 and 60 — 400 models. 
These results are similar to those found in B05 over the 
redshift range 2 < 2 < 4. On comparing the 15 — 200 and 
15 — 400 simulations which have different mass resolutions, 
with increasing resolution T_i2 is lowered by 14 per cent at 
2 = 5 and 26 percent at 2 = 6. There is an even larger drop 
of 26 per cent and 41 per cent when comparing the 15 — 100 
and 15 — 200 models at the same redshifts. In comparison, 
the drop in r_i 2 found between the 30 - 200 and 30 - 400 
models is 26 per cent at 2 = 5 and 38 percent at 2 = 6. 
Again, the trend in these results is similar to that found 
in B05, in the sense that the inferred F_i2 decreases with 
increasing mass resolution. 

Based on these results, we adopt corrections for the box 
size and mass resolution which lower the F_i2 inferred from 
the fiducial 15 — 200 simulation by 24 per cent at 2 = 5 and 
38 per cent at z = 6. Note, however, that the simulations 
are still only marginally converged. Therefore, we also adopt 
an additional numerical uncertainty of 10 per cent in our fi- 
nal analysis. From this point onwards all values of F_i2 are 
quoted including these box size and mass resolution correc- 
tions. 



2.4 The dependence of T 12 on cosmological and 
astrophysical parameters 

We now determine how F_i2 scales with various cosmolog- 
ical and astrophysical parameters at 2 = 5 and 6 using the 
simulations listed in Table 1. A full discussion of the method 
used to determine these scaling relations is outlined in B05. 
For brevity, in this work we shall simply report the results of 
this procedure. The scaling relations for T_i2 at 2 = 5 and 6 
are listed in Table 4. We set r_i 2 (2) oc A xl - Z \ where the pa- 
rameters A are listed in the top row of Table 4 and the rele- 
vant indices, x(z), are listed in the subsequent columns. Note 
that, as in B05, we have adopted F_i2 oc h 3 and T_i2 cx Vt^ 
independently of redshift. Additionally, the scaling of T_i2 
with the slope of the effective equation of state, 7, is derived 
in post-processing by pivoting the fiducial 15 — 200 model 
temperature-den sity relation around the mean gas density 
l|Viel et al.ll2004h . This will not self-consistently model any 
dynamical effects the change in temperature has on the gas 
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distribution, but it will model thermal broadening and the 
change in the neutral hydrogen fraction correctly, thus pro- 
viding a reasonable approximation. We shall use these de- 
rived scaling relations to estimate the uncertainty on T_i2 
inferred from our simulations using the fiducial parameters 
and uncertainties listed in Table 3, while also including the 
corrections for box size and resolution as discussed above. 



3 FROM IONIZATION RATE TO MEAN FREE 
PATH 

In an inhomogeneous IGM, the mean free path of an ionizing 
photon is related to the average separation between Lyman 
limit systems and the cumulative opac ity of intervening, 
lower column density H i ab sorbers ( e.g. iMeiksin fc Madaul 
1 1993! ; iMiralda-Escudel [20031 ). However, one may reasonably 
approximate the ionizing photon mean free path by the 
mean separation between Lyman limit systems alone. In- 
cluding the opacity of intervening H i absorbers will reduce 
the mean fre e path for Lyman limit photons by around a 
factor of two (|Miralda-Escude1l2003l ; IFurlanetto fc~Ohll2005l ) ■ 
but a fraction of photons at the Lyman limit will also prop- 
agate further than one mean free path. This will also be true 
for higher energy photons, reducing the effect of the inter- 
vening opacity further. Nevertheless, strictly speaking this 
assumption should be regarded as providing an upper limit 
to the mean free path. 

The mean free path model w e adopt is based on 
the po st-overlap reionization model of lMiralda-Escude et al.l 
(200(j ) (hereafter MHR00). The MHR00 model assumes the 
mean free path for hydrogen ionizing photons, A m f p , corre- 
sponds to the mean free path between regions with over- 
density A > Ai. The connection between A m f p and Ai is 
provided by 



Amfp = A (l + z) [1 - F V (A < Ai) 



-2/3 



(2) 



where Fv(A < Ai) is the fraction of the IGM volume with 
A < Ai, \qH{z) = 60 km s _1 and A m f p is expressed as a 
comoving length. MHR00 set Fv(A < Ai) using a fitting 
formulae for the volume weighted probability density func- 
tion (PDF) of the gas overdensity: 



P v (A)dA = Aexp 



(A 



-2/3 



Co) 



2(2«5 /3) 2 



A _/3 dA. 



(3) 



The fit parameters were deri ved by MHR00 at z = 2, 3 
and 4 from the simulations of Miralda- Escude et al.l {l996), 
and an extrapolation was also given at z = 6, assuming 
So = 7.61/(1 + z) and a power-law slope for the high den- 
sity tail of P — 2.5. We have checked this extrapolation at 
z = 6 against Py(A) derived from our 15-400 simulation, 
and find reasonable agreement at low A. However, even our 
high-resolution simulations underpredict the amount of gas 
at high densities due to limited resolution. Therefore, in this 
work we shall persist in using the Py(A) fits of MHR00, 
which include an analytical approximation for the high den- 
sity tail of the density PDF as a best guess at correcting for 
insufficient numerical resolution. We also extrapolate Pv(A) 
to z — 5 and z > 6 by assuming So = 7.61/(1 + z), f3 = 2.5 
and setting A and Co such that the total volume and mass 
is normalised to unity. 



In order to link the value of Ai, and hence A m f p , to a cor- 
responding value for the metagalactic hy drogen ionization 
rate, w e follow the argument presented by IFurlanetto fc Obi 
( 2005). Assuming the typical size of an absorber with over- 
density A is the local Jeans length, then in ionizati on equi- 
libriu m its column density is approximately given by (|Schavd 
l200lfl 



4.5 x 10"^ f-^r (— ) 

r_i 2 vio 4 kj V 7 J 



l + z\ 9 ' 2 



(4) 



The absorber will become optically thick to Lyman limit 
photons once TVhi > 1/cm = 1.6 X 10 17 cm -2 , where cthi 
is the hydrogen photo-ionization cross-section at the Lyman 
limit. Adopting this value for Nhi and rearranging equa- 
tion ((4]) yields the expected overdensity threshold for a Ly- 
man limit system immersed in the metagalactic radiation 
field 



Ai ~ 49.5 



T 



10 4 K 



l + z 
7 



r 



2/3 



(5) 



Therefore, given an independent constraint on F_i2 we may 
estimate Ai and hence, using the MHR00 model, the ex- 
pected mean free path for hydrogen ionizing photons propa- 
gating through the IGM. Note, however, that the mean free 
path in the IGM can only estimated in this way if it is larger 
than the mean separation between ionizing sources (i.e. once 
individual ionized regions have overlapped) . Prior to overlap 
the sizes o f the ionized regions themselves set the mean free 
path (e.g. iGnedin fc Fanll2006l ). Furthermore, before over- 
lap F_i2 will vary spatially within each ionized region, and 
also from one ionized region to the next, depending on local 
conditions. Therefore the assumption of a monotonic rela- 
tionship between F_i2 and A m f p will no longer hold. How- 
ever, at z < 6 the mean free path should still be larger or 
comparable to the typical ionizing source separation. 

The relationship between A m f p and T_i2 at various red- 
shifts is shown as the straight lines in Fig. [2] The model 
is also compared to some independent observational con- 
straints on A m f p and T_i2. The ionization rate constraints 
at z — [2, 3, 4] are taken from B05, where log Thi = 



-n.89l8:S, 



-12.07±o:i5, -12-Olto:iIl- The mean free path 



estimates are based on the fit to the observed number of Ly- 
m an limit systems per unit redshi ft at 0.40 < z < 4.69 made 
bv IStorrie-Lombardi et all (|l994T ). This yields dNhhs/dz = 
[1.5, 2.3, 3.3] at z = [2, 3, 4], where 



A 



mfp 



( 



dN L 



(6) 



H(z) \ dz J ' 

Based on fig. 2 of lStorrie-Lombardi et all l| 19941 ). we assume 
an uncertainty of ±25 per cent for the estimate of dN^s/dz 
at all redshifts. The observational data and the theoretical 
curves are consistent with each other within the uncertain- 
ties, although the observed number of Lyman limit systems 
appears to evolve slightly more rapidly. 



4 THE CONTRIBUTION TO THE IONIZING 
BACKGROUND FROM OBSERVED 
GALAXIES AND QUASARS 

We now proceed to compute the expected contribution 
to the ionizing background from LBGs and quasars at 
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Figure 2. The dependence of the comoving mean free path for 
ionizing photons on the metagalactic hydrogen ionization rate. 
The straight lines correspond to the predictions made by the 
model discussed in Section 3 at several different redshifts. In- 
dependent observational constraints on A m f p and T_i2 at z = 
[2,3,4], taken from Storrie-Lombardi et al. (1994) and Bolton et 
al. (2005) respectively, are shown by the filled symbols with error 
bars. 



2 = 5 and 6. The comoving monochromatic emissivity, e„ 
[erg s _1 Hz -1 Mpc -3 ], of a particular galaxy or quasar pop- 
ulation can be computed from the observed luminosity func- 
tion 4>(L V , z) by 



L V 4>(L V , z)dL u , 



(7) 



where in general L m i n is the luminosity corresponding to the 
limiting magnitude of the relevant galaxy or quasar survey; 
lower values require the uncertain extrapolation of the faint 
end of the luminosity function. The comoving emissivity is 
then related to the metagalactic hydrogen ionization rate by 



= A mfp (l + z) /esc / e v —^-dv. 



(8) 



assuming the comoving mean free path, A m f p , for all ioniz- 
ing photons emitted by a population of Poisson distributed 
sources is much smaller than the horizon scale. This will be 
a reasonable approximation at the redshifts we consider (see 
Table 7 for details). Again note that we assume all ionizing 
photons propagate for only one mean free path, equivalent 
to the mean separation between Lyman limit systems. Here 
/ csc is the escape fraction of ionizing photons, a v is the pho- 
toionization absorption cross-section, h p is Planck's constant 
and i^L = 13.6 eV/h p is the hydrogen Lyman limit frequency. 
Under these assumptions, the comoving emissivity is related 
to the specific intensity of the ionizing background, ft, at 
frequency v by 



ft 



A m fp 



4-k' 



(9) 



4.1 Quasar emissivity 

Firstly we consider the expected ionizing emissivity from 
quasars. To estimate the comoving quasar emissivity at 
z = 5 and 6 we have used the maximum l ikelihood fits to 
the quasar luminosity function derived by Mei ksinl (|2005T ) 
assuming a pure luminosity evolution, 



4>(L,z) 



<t>*/L«{z) 



[L/L.(*)]-A+[L/L.(*)]HV 



(10) 



where ft = -1.24, ft = -2.70, log(ft7Gpc" 3 ) = 2.51, 
log(L,(5)/L ) = 11.89 and log(L,(6)/L ) = 11.53. These 
fits are in good agreement with th e faint and bright end 
slopes repor ted bv lHunt et al.l (|2004l ) at z ~ 3 and lFan et al.l 
|200ll . 12004 ) at z > 5, respectively. However, recent work 
based on an analysis of deep X-ray data indicates the 

faint end slope at z ~ 6 may be significantly steeper 

llShankar fe Mathudl2007h . We adopt a broken power law 
for the quasar spectral energy distribution: 



oc 



-0.5 



(1050 < A < 1450 A), 



I/" 1 ' 5 (A< 1050A). 



(11) 



Assuming the photoionization cross-section has a frequency 
dependence a„ — 6.3 x 10~ la (v/i^)~ 3 cm 2 , we may then 
integrate equation ||8} to obtain: 



: 0.04 e§ 4 



where 



£l/10 24 erg s _1 Hz -1 Mpc -3 is the comoving 
quasar emissivity at the Lyman limit, a B is the spectral index 
at A < 912 A and A m f p is expressed in comoving Mpc. We 
have assumed that all ionizing photons emitted by quasars 
escape into the IGM. 

There are three main parameters on which rl 12 de- 
pends, the ionizing photon mean free path, the spectral in- 
dex and the quasar emissiv ity. As reiterated recently by 
ISrbinovskv fc Wvithd l|2007i ). the emissivity depends sen- 
sitively on the minimum quasar luminosity, L m i n , used 
when integrating equation 0. For this work we shall es- 
timate the quasar emissivity by integrating equation Q 
to a faint limit of Mab(1450) = —22, similar to the 
faintest magnitudes probed at these re dshifts by combined 
deep X-ray and optical s urveys (e.g. iBarger et all 120031 ; 
iDiikstra fc Wvithe I 120061 : IShankar fc Mathurl I2007T ). This 
corresponds to log(L m i n /£ ) = 11.17 within our adopted 
model. The comoving quasar emissivities at the Lyman 
limit computed using equation (7) at z = [5, 6] are then 
e 24 ~ [0.57, 0.21] . These em issivities are similar to those 
reported by iMeiksirj l|2005h . These values may be some- 
what larger if qua sars have harder spectra at these redshifts. 
IScott et al.l (|2004T ) find a harder spectral index of a B — —0.56 
may be appropriate below the Lyman limit, although this 
value is derived from quasars at z < 1. 



4.2 Galaxy emissivity 

We compute the expected comoving emissivity from LBGs 
in a similar fashion to the quasar emissivi ty. LBG luminos - 
ity functions are commonly expressed as a ISchechterl (|l976l ) 
function, where the number of galaxies per unit comoving 
volume is given by 
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(13) 



For the LBG luminosi t y fun ction at z = 5, we use the recent 
data of lYoshida et all feoOfJ) at (z) = 4.7, where a = -2.31, 
log(</>,/Gpc -3 ) = 5.76 and log(L*/L ) = 10.79. At z = 6, 
we consider two differing, inde pendent constra i nts o n the 
LBG luminos i ty fu nction from Bouwens et alJ (I2006T ) and 
iBunker et all (|2004h . The lBouwens et alJ (12009 ) luminosity 
function parameters are a = — 1.73, logf^/Cpc" 3 ) = 6.31 
and log(L»/L Q ) = 10.50. The lBunker et all i|2004l ) datfQare 
a = -1.6, log(0,/Gpc- 3 ) = 5.36 and log(L,/L ) = 10.75. 

We adopt a s pectral energy distrib ution based on the 
model spectrum of lLeitherer et all (|l999T ) for a galaxy of age 
500 Myr with a continuous star formation rate, a Salpeter 
IMF and metallicity Z = 0.2Z Q : 



(912 < A < 3000 A), 
(A < 912 A), 



(14) 



with an additional b reak in the spectrum at the Lyman limit, 
e L = e(1500)/6 (e.g. iMadau et al.lll999l . hereafter MHR99). 
Once again, integrating equation (0 we obtain 



: 0.03 



m~'(£r) 



A m fp \ 

40 Mpc J 



(15) 



where / osc is the escape fraction for ionizing photons, ef 5 = 
e|/10 25 erg s _1 Hz -1 Mpc -3 and A m f p is again expressed in 
comoving Mpc. 

In this instance there are four parameters which in- 
fluence estimates of F^: the mean free path, the galaxy 
emissivity, the spectral index and additionally the escape 
fraction. We attempt to take into account the uncertainty 
in the emissivity by integrating equation $7} to two differ- 
ent lower luminosity limits, roug hly corresponding to the 
faint end limits of the surveys bv lYoshida et alJ J2006) and 
iBouwens et~ai] (|2006h . These limits are M AB (1350) = -20 
and —18 respectivel y. Using these limits , we obtain e| 5 = 
0.96 and 4.31 for the lYoshida et all (|2006l ) d ata at z = 5. At 
z = 6 , we obtain ef 5 = 0.63 and 2.36 for the I Bouwens et al 
20061 ) data and ef 5 = 0.23 and 0.51 for the IBunker et al 



2004) data. These results are summarised along with the 



derived quasar emissivities in Table 5. Lastly, note that the 
spectral shape of the ionizing emission from high redshift 
galaxies is rather uncertain. For example, a harder spec- 
trum characterised by Population-Ill sta rs, a s ~ 1 (e.g. 
iBromm et al.l [200ll : iTumlinson et al.1120031 ). would increase 
the derived ionization rates by a factor of 1.5. 



5 SPATIAL FLUCTUATIONS IN THE 

IONIZING BACKGROUND APPROACHING 
REIONIZATION 

5.1 A simple model for ionizing background 
fluctuations 

At z < 4 spatial fluctuations in the metagal actic hydro- 
gen ionization rate are expected to be small (|Croft et al.l 



Fit taken from table 13 in lBouwens et at (2006) 



Table 5. Summary of the quasar and LBG Lyman limit emissiv- 
ities derived using the luminosity functions and spectral energy 
distributions given in Section 4. The last column lists the adopted 
faint end magnitude limit, M^(1350), used when integrating the 
LBG luminosity function. 



z e 24 


e s 
€ 25 


LBG data 


magnitude limit 


5 0.57 


4.31 


Yoshida et al. (2006) 


-18 




0.96 




-20 


6 0.21 


2.36 


Bouwens et al. (2006) 


-18 




0.63 




-20 


6 0.21 


0.51 


Bunker et al. (2004) 


-18 




0.23 




-20 



1 19991: iMeiksin fc White] 120041 ; ICroftl 120041 : but see also 
iMaselli fc Ferrara 1 12005). At these redshifts, the mean free 
path for ionizing photons is substantially larger than the 
mean separation between ionizing sources, therefore a spa- 
tially uniform ionizing background is expected to be a rea- 
sonable approximation. However, towards higher redshifts, 
and especially approaching the tail end of the hydrogen 
reionization epoch, spatial fluctuations in the ionizing back- 
ground are gradually amplified. This amplification is at- 
tributable to diminishing sources numbers, a smaller mean 
free path and the inhom ogeneous distribution of the ionizing 
sourc es themselves (e.g. iBolton et aUl2006l ; IWvithe fc Loebl 
I2006J) . 

However, when obtaining T_i2 from our simulations we 
assume a spatially uniform ionizing background. We must 
therefore account for the systematic error this assumption 
will impart to our measurements of r_i2. A correct numer- 
ical treatment of ionizing background fluctuations requires 
a prescription for cosmological radiative transfer, which is 
beyond the scope of this paper. We instead use a modi- 
fied version of the fl uctuating ionizing background model of 
IBolton et~ al. (2006), originally used to analyse spatial fluc- 
tuations in the quasar dominated He 11 ionizing background 
at 2 < z < 3, combined with our 30-400 simula tion. The de- 
tails o f the model implementation are given in IBolton et al.l 
(2006), to which we refer the reader. Here we only note the 
modifications made to the model for this work. 

Firstly, we now use the model to compute fluctuations 
in the ionizing background for H 1 rather than He 11 ionizing 
photons. Therefore, galaxies as well as quasars are included 
in the model by us ing the LBG l uminosity func t ion o f 
Yoshida et all i|2006h at z = 5 and IBouwens et al.l (|2006f ) 
at z = 6. The hydrogen ionization rate at some position ro 
may be written as 



AT 

Im = E 

j=i 



L 3 V fa 



47r|rj — ro| 2 h p v 



(16) 



for |rj — ro| < A m f p /(1 + z), where |rj — r () | is the proper 
distance of ionizing source j with luminosity from ro and 
N is the total number of ionizing sources brighter than Lum 
in some fixed volume. We set / csc = 1 for quasars, / csc = 0.1 
for galaxies and assume that all ionizing photons propagate 
a distance corresponding to one mean free path. The total 
number of ionizing sources brighter than absolute magnitude 
Mum, corresponding to Ln m , within a volume V at redshift 
z is then fixed by 
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Figure 3. Left: The ratio of the hydrogen ionization rate inferred from hydrodynamical simulations with a spatially fluctuating ionizing 
background to that obtained assuming a spatially uniform ionizing background at z = 5. The lower horizontal axis denotes the mean free 
path used in each fluctuating background model, while the upper horizontal axis corresponds to the ratio of the mean free path to the 
mean separation between galaxies in each model. The data are constructed from galaxy luminosity functions integrated to two different 
limiting magnitudes, Mab(1350) = —18 and —20, shown by the black and grey bars respectively. Right: As for left panel but at z = 6. 



N{z, M < Mux, 



V 



<j>{M, z) dM, 



(17) 



where <f>(M, z) is the relevant quasar or LBG luminosity 
function expressed in terms of magnitudes. The LBG lu- 
minosity functions are integrated to two different limiting 
magnitudes, M A b(1350) = -18 and -20. The brighter lim- 
iting magnitude produces fewer galaxies within a fixed vol- 
ume, increasing the amplitude of the spatial fluctuations. 
The limiting magnitude for the quasar luminosity function 
is taken to be Mab(1450) = —22 as before. 

Haloes are identified within the 30-400 simulation vol- 
ume using a friends of friends halo finding algorithm with 
a linking length of 0.2. Ionizing source luminosities are as- 
signed to the identified haloes by Monte-Carlo sampling the 
quasar and LBG luminosity functions. The haloes are pop- 
ulated by associating the brightest sources with the most 
massive haloes in the simulation volume. A more detailed ap- 
proach is impossible in this instance due to the limited reso- 
lution of the simulation, but this approach should model the 
spatial distribution of the ionizing sources reasonably well. 
The galaxy and quasar spectral energy distributions adopted 
are given Section 4 of this paper. The fluctuating ionizing 
background models are then constructed at z — 5 and 6 us- 
ing five different values for the mean free path, A m f p = 10, 
25, 50, 100 and 150 comoving Mpc. The effect of finite source 
lifetimes and light cone effects (e.o. lCroftlliool ) are not in- 
cluded in our model, but their impact on the fluctuation 
amplitude for the sm all mean free paths co nsidered here is 
expected to be small (|Wvithe fc Loebl 120061 ). 



5.2 The impact of ionizing background 

fluctuations on T_i2 inferred from simulations 

The metagalactic ionization rate we infer from our simu- 
lations using a spatially fluctuating ionizing background, 
rLi2i i s compared to F_i2 obtained assuming a spatially 
uniform ionizing background in Fig. [3] The left panel shows 
the 2 = 5 data, with r'_ 12 /r_i2 for each fluctuating ioniz- 
ing background model plotted as the black bars for a lim- 
iting magnitude of Mab(1350) = —18 and grey bars for 
Mab(1350) = —20. The lower horizontal axis corresponds 
to the comoving ionizing photon mean free path, while the 
upper horizontal axis shows the ratio of the mean free path 
to the mean separation between galaxies in the model, (Z g ). 
Similar data at z = 6 are shown in the right hand panel of 
Fig. [3] Note that (l s ) is an average galaxy separation com- 
puted assuming the galaxies are uniformly distributed. In 
reality the inhomogeneous distribution of ionizing sources 
within our model results in underdense regions being sur- 
rounded by fewer ionizing sources than overdense regions, 
further increasing the amplitude of the spatial fluctuations. 

As one might expect, the impact of spatial fluctua- 
tions is largest when A m f p /(Z g ) is small. Increasing the num- 
ber of ionizing sources within one mean free path of a 
fixed point (either by increasing A m f p or adopting a fainter 
limiting magnitude) reduces the amplitude of the fluctua- 
tions and thus the departure of F'_i2 from the ionization 
rate inferred assuming a spatially un iform ionizing back- 
ground. As found in previous studies ( Gncdin fc Hamiltonl 
|2002| ; iMeiksin fc White! |2004| ; iBolton et al.l l2006h the as- 
sumption of a spatially uniform ionizing background tends 
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to underestimate the ionization rate required to reproduce 
the observed Lya effective optical depth at high redshifts. 
The presence of large ionization rates in regions where the 
Lya opacity of the IGM is already rather low or very high 
increases the mean F_i2 but has little impact on the ob- 
served effective optical depth. For plausible values for the 
mean free path at z — 5 this effect is around the 5 — 10 per 
cent level. 

Interestingly, there is a different trend at z = 6 for 
progressively smaller values of A m f p . The difference between 
rli2 and T_i2 is reduced as A m f p decreases, and T_i2 is 
actually larger than r'_ 12 for the model with the smallest 
value of Amfp/^g) at z = 6. This result was also noted by 
iMeiksin fc White! i|2004h . In this instance, the volume of the 
IGM where there are no ionizing sources within a mean free 
path become more common, reducing the spatially averaged 
ionization rate. 

In summary, we find the ionization rate inferred from 
simulations of the IGM which assume a uniform ionizing 
background at z = 5 and 6 is usually underestimated. How- 
ever, as regions which are isolated from ionizing photons 
become more common, the T_i2 can actually be overesti- 
mated when assuming a uniform ionizing background. For 
plausible values of A m f p , spatial fluctuations in the ioniz- 
ing background induce a systematic offset of about +10 per 
cent and an uncertainty of about ±10 per cent to the T_i2 
we infer from our simulations at z — 5 and 6, respectively. 
We shall include these uncertainties when quoting our final 
estimate of T_i2 in Section [6] 



6 THE METAGALACTIC HYDROGEN 
IONIZATION RATE AT Z=5 AND 6 

6.1 Constraints on T_i2 from the IGM 
Lya opacity 

We summarise our constraints on the metagalactic hydro- 
gen ionization rate at z — 5 and 6 in Fig. [4] The values of 
r_i2 which reproduce the IGM Lya effective optical depth 
are plotted as the filled triangles with la error bars. Our 
central estimates and their uncertainties are based on T_i2 
inferred from our fiducial simulation (15-200), rescaled to 
correspond to the fiducial parameter values and ranges listed 
in Table 3. We do this using the scaling relations we derived 
for r_i2 in Section 2. The total uncertainty on F_i2 is then 
computed by adding the contributions from individual pa- 
rameter uncertainties in quadrature. The total error budget 
is summarised in Table 6. We find the Lya effective optical 
depth is reproduced by logFin = — 12.28^0 23 at z = 5 and 
logFrn < -12.72 at z = 6. 

These data are compared to other observational con- 
straints taken from the literature in Fig. [4] Th e open cir- 
cles w ith error bars correspond to the data of iFan et al71 
(2006), derived from the IGM Lya effective optical depth 
at z ^5 — 6. Rather than using hydrodynamical simula- 
tions, ^^et^aT] (|2006l ) compute T_i2 u sing the fluctuatin, 
Gunn-Peterson appr oximation (FGPA. iRauch et al.l K199 
IWeinberg et al.lll999T ) combined with the IGM density dis 
tribution derived by MHR00. The data are in agreement 
with our constraints at z = 6, although our upper limit 
on r_i2 is somewhat higher. However, at z = 5 our de- 
termination of F_i2 is significantly larger than the data of 



Table 6. The error budget for T_i2 based on the fiducial val- 
ues for various cosmological and astrophysical parameters given 
in Table 3. The uncertainties due to the impact of ionizing back- 
ground fluctuations and marginal numerical convergence are also 
listed. The total error is obtained by adding the individual errors 
in quadrature. 
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z = 5 (per cent) 


z = 6 (per cent) 


To 
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+44/-32 


r2 m 


+24/-17 


+26/-18 
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+24/-19 


o"8 
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Fluctuations 


+10 


±10 


Numerical 


±10 


±10 


Q h h 2 
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±6 


±6 
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+66/-40 
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Figure 4. Evolution of the metagalactic hydrogen ionization rate 
over the redshift range 2 < z < 6. Our new estimates of r_i2 are 
shown as the filled triangles at z = 5 and 6. Other data points 
based on the IGM Lya effective optical depth are taken from 
Tytler et al. (2004) (star), Bolton et al. (2005) (open diamonds) 
and Fan et al. (2006) (open circles). The box corresponds to the 
constraints on T_i2 obtained from the proximity effect by Scott 
et al. (2000), and the filled diamond is the recent estimate of 
r_i2 based on emission from Lyman break galaxies at z = 3 
(Shapley et al. 2006). The solid and dotted lines correspond to the 
ionizing background models of HM01 for, respectively, galaxies 
and quasars and quasars only. 



IFan et al. I (|2006l '). This is not unexpected; the FGPA should 
be used with caution, as in practice the dependence of the 
Lya optical depth (not to be confused with r e g) on vari- 
ous parameters is a ctually rather different to the cano nical 
FGPA scaling (B05. lTvtler et al.ll2004l : ljena et al.ll2005h . We 
therefore expect our constraints on F_i2, based on our ex- 
tensive suite of hydrodynamical simulations, to be more ro- 
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bust. The curves in Fig. [4] correspond to the F_i2 predicted 
by the ionizing background models of HM01 for galaxies and 
quasars (solid line) and quasars only (dotted line). These 
models are described fully in B05. Our data are in good 
agreement with the model based on the combined emission 
from galaxies and quasars, which is consistent with the inter- 
pretation that quasars alone cannot maintain the IGM in its 
highly ionized state at z = 6 |Miralda-Escude 2003; Meiksinl 
120051 : ISrbinovskv fc Wvithdl2007l ). We also compare these 
curves to lower redshift constraints on T_i2 based on the 
Lya effective optical depth. The open d iamonds are the dat a 
from B05 and the star is taken from iTvtler et alJ (|2004h . 
Measurements of r_i2 can also be made using the quasar 
proximity effect (e.g. 



Baitlik et ail 1 19881 ; iGiallongo et all 



1 19961 : IScott et aD 120001 ) , although these estimates can be 
subject to a bias int roduced by the overdense regions in 
which quasars reside jRollinde et al.ll2005l ; iGuimaraes et all 
120071 : iFaucher-Gieuere et al.ll2007l ) ~The box in Fig. [4] cor- 
responds to the r_i2 inferred by IScott et all (|200oF All 
data are consistent with the ionizing background having a 
substantial contribution from young star forming galaxies. 
T he filled diamond als o shows the recent direct estimate 
bv lShaplev et all (|2006l ) of the contribution of LBGs to the 
ionizing background at z ~ 3. The measurements are con- 
sistent with around 50 per cent of the ionizing photons at 
z ~ 3 originating from galaxies. 

Lastly, we compare our data with previous constraints 
on r_i2 obtained from other simulatio ns of the Lya for- 
est. |McDonaid^Miral^Escudi i|200l[ ) use the hydro-PM 
code of iGned in fc Huil p998) to run a simulation with a 
box size of 8.9/i _1 comoving Mpc and a 256 3 grid. They 
obtain T_i2 = 0.13 ± 0.03 at z = 4.93 assuming a flat 
ACDM cosmological model with fi m = 0.4, Q. h h 2 = 0.02, 
h = 0.65, a 8 = 0.79 and n = 0.95 for r cff = 2.3, T = 
2 x 10 4 K and 7 = 1. We rescale their constraint on T_i2 to 
match our fiducial cosmological parameters, listed in Table 
3, with the scaling relations derive d at z = 5 in Table 4. 
iMcDonald fc Miralda-Escudel l|200ll) do not apply a numer- 
ical resolution correction factor to their data, so for a fair 
comparison with our data we also multiply their rescaled 
r_i2 constraint by the appropriate resolution correction fac- 
tor of 0.76 we obt ained at z = 5. This p rocedure yields 
r_i2 = 0.44 ± 0.10. iMeiksin fc White! (120041 ) employ a simi- 
lar pseudo-hydrodynamical approach with a simulation box 
size of 25/i -1 comoving Mpc, 512 3 dark matter particles and 



a 1024' 1 force mesh. They obtain r_i 2 = 0.3F 



at z = 5 



for a flat ACDM model with Q m = 0.3, VL h hT = 0.022, 
h = 0.70, as = 0.92, n = 0.95 and r cff = 2.12. The gas 
is assumed to follow an effective equation of state with 
To = 2 x 10 4 K and 7 = 1.5. Rescaling their result yields 
r_i2 = 0.52^q']^. Both of these results are good agree- 
ment with our constrain t of T_i2 = 0.52 j^o |i ■ At z = 6, 
ICen fc McDonald! l|2002l) use an Eulerian hydrodynamical 
simulation with a box size of 25/i -1 comoving Mpc and 
a 768 3 grid to obtain a la upper limit of T_i2 < 0.096. 
They adopt fi m = 0.3, Q, h h 2 = 0.02, h = 0.67, a 8 = 0.9 
and n = 1 with r e a = 5.57; no values for To or 7 were 
quoted. Rescaling their constraint on r_i2 as before, and 
including an appropriate resolution correction factor of 0.62 
at z = 6 yields T_i2 < 0-11 a ssuming To = 10 4 K and 
7 = 1.3. IMeiksin fc White! (|2004T ) quote a la upper limit of 
r_i2 < 0.14 for Tefi > 5.12 at z = 6 . Rescaling their result 



gives F_i2 < 0.20. These results are again consistent with 
our upper limit of F_i2 < 0.19. 

6.2 Constraints on T_i2 from the observed galaxy 
and quasar populations 

We now compare our estimates for F_i2 at z = 5 and 6 to 
the ionization rates estimated from the Lyman limit emis- 
sivities listed in Table 5. The total ionization rate is F_i2 = 
T^ 12 +ri 12 , where the quasar and galaxy contributions are 
computed using equations (I12|l and (|15|) . This is plotted as 
a function of A m f p in Fig. [5] for four different values of the 
LBG escape fraction, / csc . The escape fraction for ionizing 
photons is very uncertain, a lthough recent work suggests 
/esc may vary with redshift ( Ra zoumov fc Sommer-Larsenl 
2006), and that / esc > 0.1 may be appropriate at z ~ 6 
jlnoue et afll20odh 



The top left and right pane ls in Fig. [5] d isplay 
F_i2 obtained by integrating the lYoshida et al.l (|2006l ) 
LBG luminosity function to Mab(1350) = —18 and 
M ab (1350) = —20. resp ectively. The limiting magnitude of 
the lYoshida et al.l (120061 ') data is roughly M AB (1350) = -20. 
The quasar contributions to F_i2 correspond to 15 per cent 
for M AB (1350) = -18 and 44 per cent for M AB (1350) = 
—20, assuming / csc = 0.1. For comparison, the thick solid 
lines show our determination of T_i2 from the IGM Lya ef- 
fective optical depth at z = 5, with our corresponding 
estimate for A m f p using the model described in Section 
3. For r_i2 = 0.52^0 21 we obtain a mean free path of 
Amfp = 114^34 comoving Mpc. The combined emission from 
galaxies and quasars at z — 5 is consistent with the IGM 
Lya effective optical depth, assuming Mab (1350) = —18 
and / csc > 0.15, or M A b(1350) = -20 and / csc > 0.6. 

The central panels show similar data for r_i2 at z = 6 
using the LBG luminosity function of lBouwens et al.ll|2006h . 
A constant quasar contribution is again added, correspond- 
ing to 11 percent of the total T_i2 assuming a magni- 
tude limit of Mab(1350) = —18 and 31 per cent for 
M AB (1350) = -20, both for / e 3C = .1. The limiting mag- 
nitude of the iBouwens et al.l l|2006h data roughly corre- 
sponds to Mab(1350) = —18. The thick solid line gives 
our upper limit of F_i2 < 0.19 and A m f p < 37 comoving 
Mpc. Galaxies and quasars provide the required number of 
photons needed to reproduce the IGM Lya effective opti- 
cal depth for M A b(1 350) = -18 and / csc > 0.2. For the 
IBouwens et al. (2006) data there does not appear to be any 
deficit in the number of ionizing photons required to ionize 
the IGM at z = 6. A similar conclusion was also reached by 
IBouwens et~ai1 (|2006l ). 

The bottom panels show T_i2 obt a ined from the LBG 
luminosity function of iBunker et all (|2004h at z = 6. 
A constant quasar contribution of 36 percent of the to- 
tal r_i2 is added for Mab(1350) = —18 and 55 per 
cent for Mab (1350 ) = -20, assuming / csc = 0.1. The 
IBunker et al.l (|2004l ) luminosity function has a substantially 
lower galaxy number densi ty at the faint end compared 
to the IBouwens et al] (120061 ) data. Consequently, the T_i2 
computed from the LBG and quasar emissivities require 
a much higher value of / OS c for consistency with our con- 
straints on r_i2 and A m f p from the IGM Lya effective op- 
tical depth. At Mab (1350) = —18, roughly corre s pond - 
ing to the limiting magnitude of the IBunker et al.l (|2004f ) 
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Figure 5. Comparison of our constraints on T_i2 and A m f p obtained from the IGM Ly« effective optical depth at z = 5 and 6 (thick 
solid lines) to T_i2 computed from the Lyman limit emissivity of the observed galaxy and quasar populations at the same redshifts (thin 
straight lines). These are plotted as a function of A m f p for four different values of the ionizing photon escape fraction from LBGs. Top 
panels: Data at z = 5 computed using the LBG luminosity function of Yoshida et al. (2006). The left and right hand panel correspond to 
the LBG emissivities computed by integrating the luminosity function to limiting magnitudes of Mab(1350) = — 18 and —20 respectively. 
Central panels: Data at z = 6 computed from the LBG luminosity function of Bouwens et al. (2006). The left and right hand panels 
correspond to the galaxy emissivities computed by integrating the luminosity function to limiting magnitudes of Mab(1350) = — 18 and 
—20 respectively. Bottom panels: As for central panels but using the galaxy luminosity function measured by Bunker et al. (2004). 
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galaxy luminosity function, / csc ~ 1 is required, whereas for 
Mab(1350) = —20 LBGs and quasars cannot account for 
the required ionization rate at z — 6. An escape fraction of 
/csc ~ 1 is probably implausibly high but given the rather 
uncertain spectral energy distribution of the sources this is 
probably not a cause for concern. 

However, this time our conclusion is ra ther different 
from the one reached by iBunker et al, (2004). On the ba- 
sis of their luminosity function they concluded that there is 
a significant deficit in the number ionizing photons at z — 6 
required to keep the IGM highly ionized. Even for / csc = 1, 
Bunker et al. found the number of photons required to ionize 
the IGM provided by galaxies suffers a substantial shortfall. 
Note that this difference is not largely due to our assump- 
tion that quasars also contribute to the ionizing emissivity. 
For the assumed emissivities given in Table 5, quasars only 
increase the total number of ionizing photons by around a 
factor of two assuming an LBG escape fraction of / esc = 0.1, 
and by much smaller factors for larger escape fractions. As 
discussed in the next section, this discrepancy can be at- 
tributed to the assumption of an unduly large H n clumping 
factor when estimating the number of photons required to 
keep the IGM highly ionized. 

Lastly, one may have some concern over the above con- 
clusions regarding the adopted mean free path. Recent cos- 
mological radiative transfer simulations favour a value for 
Amfp which is lower by a around a factor of two at z — 5 
l|Gnedin|[200l ; iGnedin fc Far]|2006l ). However, halving our 
estimate for A m f p at z — 5, we still find that the observed 
galaxy and quasar population should provide enough ion- 
izing photons when assuming Mab(1350) — —18 with a 
slightly larger escape fraction of / CS c > 0.3. At z = 6, our 
A m fp is comparable to the recent simul ation data, as we ll 
as the constraints on A m f p obtained by iFan et ah ] (|2006t ). 
However, if A m f p is nevertheless reduced by a factor of two, 
one instead re q uires / eS c > 0.4 for consistency with the 
iBouwens et al.l (|2006l ) luminosity function, integrated to a 
limiting magnitude of Mab(1350) = —18. These escape frac- 
tions may be slightly high, but we think that given the 
galaxy and quasar emissivities we derive are probably lower 
limits and that the spectral shape of the sources is rather 
uncertain they are not a cause for concern. Note further 
that the re-emission of ionizing photons by recombining ions, 
which has been neglected here, may also raise the ionization 
rate by as much as 50 per cent if the spectrum of ioniz- 
ing photons is hard en ough to significantly ionize helium 
(|Haardt fe Madaull 19961 ). 



6.3 The clumping factor and the minimum 
ionizing emissivity required to balance 
recombinations 

The emission rate of ionizing photons per unit comoving vol- 
ume, -/Vrec required to balance recombinations in the ionized 
part of the IGM is given by (e.g. MHR99) 

7V r cc = ^4 = 10 5 °-°Chii (i^£) 3 s" 1 Mpc" 3 , (18) 

where (tih) is the mean comoving hydrogen density in the 
Universe and (t re c) is the volume averaged recombination 
time for ionized hydrogen with an effective H n clumping 
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Figure 6. Comparison of our constraints on T_i2 and A m f p ob- 
tained from the IGM Lya effective optical depth at z = 6 (thick 
solid line) to the minimum hydrogen ionization rate required to 
ionize the IGM up to some density threshold corresponding to 
the clumping factor Cjjii at z = 6 (thin straight lines). 



factor Chii = ("■Hii)/( n Hii) 2 - The clumping factor param- 
eterises the inhomogeneity of the ionized hydrogen in the 
IGM and therefore it is a redshift dependent quantity which 
increases with decreasing redshift. It does not include the 
contribution from collapsed, virialised objects or clumps 
that are self-shielded from the ionizing background. There- 
fore, Chii can be significantly smaller than the total bary- 
onic clumping factor. MHR99 adopt Chii = 30, based 
on the value at z = 5 compu ted from the cosm o logica l 
radiative transfer simulation of IGnedin fc Ostrikeri l|l997f ). 
IBunker et al.l (|2004h also adopt Chii = 30 when comput- 
ing the number of ionizing photons needed to keep the IGM 
highly ionized at z = 6, as do a number of other authors. 
However, a wide range of ot her values for Chii have been as- 
sumed in the literature (see ISrbinovskv fc Wvith"ell2007l for 
details). Note that we use the case-B recombination coef- 
ficient evaluated at 10 4 K in equation JT8| {e.g. MHR99). 
However, if the case- A recombination coefficient is the more 
appropriate choice, iV r cc would be raised by a factor of 1.6. 
Alternatively, since the recombination coefficient scales as 
T -0 ' 7 , a gas temperature of 2 x 10 4 K would lower 7V rC c by 
around the same factor. Lastly, we note equation (|18p can 
only predict the number of ionizing photons required to bal- 
ance recombinations when the average recombination time 
is smaller than the age of the Universe. 

The emission rate of ionizing photons can be converted 
to a flux and hence a hydrogen ionization rate by assuming 
a photon mean free path and a typical source spectral index, 
a s . The corresponding specific intensity at the Lyman limit 
is then (e.g. MHR99) 

jrcc _ —\ mlp N tec h pCts (l + zf, (19) 
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Figure 7. Observational constraints on the emission rate of ionizing photons per comoving Mpc, iVjoaj as a function of redshift. All results 
are computed assuming a s = a\, = 3 and, in the case of the LBG and Lya emitter emissivity estimates only, / eS c = 0.2. The scale on the 
right-hand vertical axis corresponds to the number of ionizing photons emitted per hydrogen atom over the Hubble time at z = 6. The 
filled triangles give an estimate of iVi Dn based on the constraints for T_i2 and A m f p obtained from the Lya effective optical depth in this 
work and in B05. The inverted triangle at z = 5 and the diamond and star at z = 6 correspond to estimates of Af; on based on the Lyman 
limit emissivities of LBGs and quasars listed in Table 5. The data have been slightly offset from their actual redshifts for clarity. An escape 
fraction of / csc = 0.2 has been assumed in this instance. At z > 6, the open squares and circles are derived from the upper limits on the 
comoving star formation rate per unit volume inferred by Bouwens et al. (2005) and Richard et al. (2006), respectively. The cross is derived 
from the number density of Lya emitters estimated by Stark et al. (2007). Three simple models for the evolution of JV; on are also shown 
as the dotted, short dashed and dot-dashed lines. The solid lines correspond to the emission rate of ionizing photons per unit comoving 
volume, N T ec, needed to keep the IGM ionized for various H II clumping factors. At 2 < z < 6, we find Ni on is characterised by a power 
law, Nicm(z) = 10 50 - 5- 0.06(z— 6) g — 1 Mpc -3 , shown by the long dashed line. Adopting a mean free path which is a factor of two smaller or 
setting a B = a^, = 1 will double the emissivity derived from the Lya forest opacity. 



Assuming the spectral index of the ionizing background, 
Qb = Q s = 3, the minimum ionization rate required to bal- 
ance recombinations in the IGM with a clumping factor Chii 



prcc 
L -12 



: 0.07 C H 



(t)( 



a h + 3 



A 



m f p 



40 Mpc 



(20) 



although note that reprocessing of the intrinsic source spec- 
trum by the IGM can alter the spectral shape of ionizing 
background somewhat l|Haardt fc MadavJ Il996l ). In Fig. [6] 
we compare this ionization rate required to balance recom- 
binations for different H n clumping factors at z — 6 (thin 
straight lines) to our determination of T_i2 and A m f p from 
the IGM Lya effective optical depth (thick solid line). Our 
result for the photoionization rate is consistent with the 
value needed to balance hydrogen recombinations at z = 6 
if Chii < 3. Adopting Chii = 30, a value of T_i2 substan- 
tially higher than that inferred from the Lya effective op- 
tical depth would be required to keep the IGM highly ion- 
ized at z = 6, as was already noted by MHR99. Even us- 
ing a value for A m f p which is a factor of two smaller, our 
constraints on T_i2 still require that Chii < 7 at z = 6. 
Comparing to some re c ent co smological radiative transfer 
simul ations, llliev et al.l ll2006t ) predict Chii < 2 at z > 11, 
while ISokasian et al.l (|2003h find Chii ~ 4 at z = 3. How- 



ever, simulations with higher spatial resolution may favour 
larger values for Chii- Nevertheless, we conclude adopting 
Chii = 30 overestimates the H n clumping factor of the IGM 
at 2 = 6 by at least factor of 4 and possibly by a factor of 
10. We therefore emphasise that care should be taken when 
drawing conclusions from equation (|18[1 as Chii is rather 
uncertain. 



7 THE IONIZING EMISSIVITY - EVIDENCE 
FOR A PHOTON-STARVED AND 
EXTENDED PERIOD OF REIONIZATION 

7.1 The ionizing emissivity at 2 = 2 — 6 

In the last section we focussed on the hydrogen ionization 
rate and the question of whether or not the expected ioniz- 
ing flux from observed high-redshift galaxies and quasars is 
sufficient to keep the IGM highly ionized at z — 5 — 6. We 
now turn again to the comoving ionizing emissivity in order 
to discuss whether or not the inferred emissivity, if extrapo- 
lated to higher redshift, is sufficient to reionize hydrogen in 
the first place. 

The filled triangles in Fig. 7 show the result of turning 
the constraints on T_i2 from the Lya opacity at z — 2 — 6 
derived here and in B05 into an emission rate of ionizing 
photons per unit comoving volume. For this we use the re- 
lation 
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(?) m 



A 



mfp 



40 Mpc 



Mpc" 



(21) 



which assumes the mean free path is much smaller than the 
horizon. Note that this expression is independent of / eac , but 
does depend inversely on the mean free path for ionizing 
photons. We have estimated A m f p at z = 2 — 6 by using the 
simple model discussed in Section 3 combined with the T_i2 
constraints from this work and B05. The assumed values are 
summarised in Table 7 with a comparison to the horizon 
scale, and these are used to calculate the filled triangles in 
Fig. 7. The error estimates on Mon computed using equa- 
tion (|21[) include the error estimates for T_i2 and the cor- 
responding error in our estimate of the mean free path. The 
comoving ionizing emissivity inferred from the Lya opac- 
ity is nearly constant over the redshift range 2 < z < 6. 
Assuming a spectral index a s = «b = 3 above the Lyman 
limit, it is well characterised by a slowly rising power law, 

N ion (z) = 10 50.5-0.06( Z -6) g -l Mpc -3 i shown by the long 

dashed line. A harder spectral index of a s = «b = 1 would 
double 7Vi on , while a smaller mean free path would increase 
it linearly. 

For comparison the inverted triangle, diamond and 
star at z < 6 show the estimates of iVi on based on the 
quasar and LBG emissivities listed in Table 5 integrated 
to Mab(1350) = —18, assuming / csc = 0.2 and a s = 3. 
As noted earlier, / csc > 0.2 is required for consistency with 
our measurements of T_i2 from the Lya effective optical 
depth and estimates of A m f p . The solid curves show jV r cc as 
given by equation (|18|l for four different time independent 
H ii clumping factors. Note the curves are cut off at the red- 
shift where the recombination time exceeds the age of the 
Universe. We have not attempted to estimate errors for these 
estimates of the ionizing emissivity, but note again that the 
assumed LBG spectral shape is rather uncertain and de- 
pends strongly on the IMF, metallicity and star formation 
history. Together with the uncertainty on the escape fraction 
this leads to an error which may well approach an order of 
magnitude. In contrast, the estimates from the Lya opacity 
are significantly more robust. The biggest uncertainties as- 
sociated with these data are the Lyman limit spectral index 
and the mean free path of ionizing photons, w hich introduce 
an uncertainty of around a factor of 2-3 (see lMeiksinll2005l 
for a detailed discussion of this point). 

Finally, it is illustrative to express M on in terms of the 
number of io nizing photons emitted per hydrogen atom per 
time interval IjMiralda-Escude 2003). This quantity is shown 
on the right vertical axis in Fig. 7 for a time period corre- 
sponding to the Hubble time at z = 6. For a s — «b = 3 the 
comoving ionizing emissivity inferred from the Lya opacity 
corresponds to the emission of ~ 1.5 photons per hydro- 
gen atom in a period that corresponds to the age of the 
Universe at z = 6. For a s — ab = l this number would 
be twice that. Reionization must therefore have occurred in 
a photon- starved manner unless the ionizing emissivity was 
much higher at z > 6. 



7.2 Three simple models for the ionizing 
emissivity of galaxies at z > 6 

We now proceed to discuss plausible extrapolations of 
the ionizing emissivity towards higher redshift. For this 
we firstly consider the recent tentative upper limits on 
the UV emissivity from searches for high-redshift galax- 
ies at z > 6. Estimates based on the Hubble Ultra Deep 
Field yield e(1500) < 4.0 x 10 25 erg s" 1 Hz" 1 Mpc" 3 at 
z ~ 7.5 and e(1500) < 1.0 x IP 25 erg s _1 Hz -1 Mpc" 3 
at z ~ 10 l|Bouwens et al.l 120051). Similar re s ults h ave 
also been recently reported by Mannucci et al.l (|2007l ) at 
2 ~ 7. Alternative constraints from near infra-red observa- 
tions around lensing clusters give higher values, e(1500) < 
Mpc" 3 at z ~ 8 and e(1500) < 
z ~ 9 (Ric hard et al.l 
2006). These data omit a correction for dust extinction and 
are computed by integrating the LBG luminosity function 
to a lower luminosity limit of 0.3L„(,z = 3). In addition, re- 
cent searches for faint gravitationally lensed Lya emitters 
at 8.5 < z < 10.2 behind for eground galaxy clusters, pre- 
sented by IStark et al.l j2007b) , have given an upper limffl 
of around £L yo .<2 x 10 41 erg s _1 Mpc -3 , assuming all six 
of their candidates are real. 

The upper limits for e(1500) can be related to the emis- 
sion rate of ionizing photons per unit comoving volume, 
iVion, by 



5.3 x 10 2e erg s" 1 Hz" 1 
7.9 x 10 26 erg s" 1 Hz" 1 Mpc -3 at 



£ L/esc _ ^q49.' 
hpCt s 



^(frObr) s " Mpc " 3 ' (22) 



where we again assume el = e(1500)/6 and a source spectral 
index of a s = 3. To convert the Lya emissivity to Mon we 
use, 



3 /, 



csc c Lyo 



2 1 ,/esc hpVj^yo 



1 r>51-3 
10 £41 



fcsc \ ( 

"oir J v 



l - u. 

0.9 



Mpc" 



(23) 



where £41 = eLya/10 41 erg s _1 Mpc -3 , fcZ a is the escape 
fraction of Lya photons and we assume case B recombina- 
tion, corresponding to the pro duction of two Lya photons for 
every three ionizing photons (|Osterbrock Ill98l)l ). We adopt 
/esc" = 0. 2 based on the l ower l imit inferred from the obser- 
vations of lGawiser et al. I (|2006h at z ~ 3. 

Three simple models for the evolution of JVi on at z > 6 
are displayed along with these data points in Fig. 7. The 
models are chosen simply to bracket a range of plausible 
evolutionary histories of iVion at z > 6, which are an- 
chored at our measurement at z = 6. Model 1, shown 
by the dotted line, is a dec lining power law consistent 
with the lBouwens et all (|2005l ) data. Model 2, correspond- 
ing to the short dashed line, assumes Ni on remains con- 
stant at z > 6. Lastly, the dot-dashed line corresponds 
to model 3, a double exponential which peaks at z = 
9. Note that recent estimates of the mass assembled in 
star forming gala xies at z ~ 6 indicate such an increase 
may be plausible l|Mobasher et al.l 120051 ; lEvles et all 120061 ; 



2 Estimated from the upper limit on the number of sources with 
a Lya luminosity brighter than L, presented in fig. 11 of Stark et 



al. (2007b). 
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Table 7. Summary of r_ 12 inferred from the Lyc? forest effective 
optical depth at 2 < z < 6. The data at z = 2 — 4 are taken 
from B05, while the 2 = 5 — 6 values are based on this work. 
The corresponding estimates for the mean free path in units of 
comoving Mpc are obtained using the model discussed in Section 
3. For comparison, we also list the horizon scale in comoving 
Mpc, r nor , computed to three significant figures for our fiducial 
cosmology. 
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IStark et all l2007al ; lYan et all 120061 ). Alternatively a dis 
tinct population of high redshift ionizin g sources such as 



2003 




2005 





population-Ill stars or mini-quasars (e.g. Venkatesa n et al 
"" Sokasian et al.ll2004l; iMadau et al.ll2004l; iRicotti et al 



Choudhurv fc Ferrarall2005l ; lMeiksinll2005l ) may be 



sponsible for an increased ionizing emissivity at high red- 
shift. For all three models we assume Mon = at z > 15 



and N ion (z) 



10 



50.5-0.06(z-6) 



s" 1 Mpc" -3 at z < 



We can now obtain the filling factor of H n in the 
IGM f or the three models for Nj on by solving (e.g . 
MHR99 lWvithe fc Loebll2003l ; IChoudhurv fc Ferrarall2005n 

dQmi Ni 



tit 



<n H > 



— QhiiChii- — rr^aB(r), 



(24) 



where cib(T) is the case-B recombination coefficient, which 
we evaluate at T = 10 4 K, {««) is the mean comoving hydro- 
gen density and a = (l + jz)~ 1 is the cosmological expansion 
factor. We take Qhii = 1 as a proxy for the end point of 
hydrogen reionization. This occurs at progressively earlier 
redshifts for either a larger value of M on or a smaller value 
of Chii- The H n filling factor is shown as a function of red- 
shift for our three simple models for the evolution of Mon in 
Fig. [8] Assuming a time independent H n clumping factor 
of Chii = 2, Qhii = 1 is reached at z = 4.16 for model 1, 
z — 4.69 for model 2 and z — 6.22 for model 3. Hydrogen 
reionization is only complete by z = 6 in model 3, and is 
the only model consistent with the obser ved flux distribu- 
tion in the spectra of high redshift quasars (|Fan et al. ~l l2006l ; 
iBecker et aLll2007f ). 

However, before proceeding we should discuss the two 
main uncertainties involved in calculating Qhii from our 
three models; the source spectral index a s and the clumping 
factor Chii- The filling factors so far were obtained assum- 
ing a sou rce spect r al ind ex of a a = 3. As discussed in some 
detail by iMeiksinl (2005), a harder spectral index will in- 
crease JVi on . As an example, adopting a s = a? D = 1 doubles 
the Mon we infer at our anchor point of F_i2 at z = 6. The 
solid line labelled model 2b in Fig. [8] shows Qhii computed 
for a constant iVj on twice that in model 2. The resulting Qhii 
evolution is very similar to model 3 and is now consistent 
with hydrogen reionization being complete by z — 6. On the 



Figure 8. The H II filling factor as a function of redshift com- 
puted using the three models for the redshift evolution of Mon 
shown in Fig. 7. A time independent H II clumping factor of 
Chii = 2 and a source spectral index of a B = 3 has been assumed 
in this instance. The solid line shows the filling factor expected 
for model 2 if Ni on is raised by a factor of two, corresponding to 
a source spectral index of a B = 1. 



other hand, adopting a softer spectral index of a B — 5 (e.g. 
iBarkana fc Loebll200ll) would push reionization to even later 
times, requiring an even larger increase in Mon to achieve 
Qhii = 1 by z = 6. Our adopted value for the clumping fac- 
tor at is less of a concern. Even if we assume a uniform IGM 
(Chii = 1), model 2 still only predicts Qhii = 1 by z = 5.26. 
In addition, the clumping factor will become larger towards 
lower redshifts. 

Finally, we note that a particularly hard spectral in- 
dex has implications for the He 11 reionization history. The 
much larger ionization threshold of singly ionized helium 
relative to neutral hydrogen ensures that He 11 reionization 
is postponed until sources with sufficiently hard spectra 
become numerous. Quasars are likely to be the primary 
source of these energetic photons, and there is some evi- 
dence that the tail-end of He 11 reionization does roughly 
coincide w ith the peak in the quasar number de nsity around 
2 ~ 2 - 3 (|Shull et alj|2004 iBolton et al.ll2006l ). However, a 
spectral index of a s — 1 at z > 6 would result in an earlier 
He 11 reionization epoch, at odds with the current observa- 
tional evidence. It may therefore be possible to rule this 
scenario out. 



7.3 The IGM neutral hydrogen fraction and 
electron scattering optical depth 

The H 11 filling factors calculated in the last section may be 
used to estimate the volume weighted neutral fraction in the 
IGM for the proposed Ni on models by solving 



(/hi)v = 1 + Qhii 



f AlGM dA/ H i(A)JV(A) 

f AlGM dA p v(A) 



(25) 
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Figure 9. Left: Evolution of the volume weighted neutral hydrogen fraction in the IGM as a function of redshift. The filled triangles and 
open diamonds at z < 6 correspond to the (/hi)v computed using the constraints on T_i2 obtained in this work and in B05, respectively. 
The solid line gives the (/hi)v calculated using the UV background model of HM01 for galaxies and quasars. Other observational constraints 
at z > 6 come from Mesinger &; Haiman (2007), Totani et al. (2006) and Kashikawa et al. (2006), shown by the inverted triangle, cross 
and square respectively. The expected evolution in (/hi)v at z > 5 for the three models for -/V; on shown in Fig. [7] are shown as the dotted, 
dashed and dot-dashed lines. Right: The cumulative electron scattering optical depth as a function of redshift computed from the three 
models for 7V; on shown in Fig. [7] 



when Qhii < 1 and 

ff IGM dA f H i(A)iV(A) 

when Qhii > 1- Here /hi(A) is the residual neutral hydrogen 
fraction in an ionized region with overdensity A and Aigm 
is the upper limit for the gas overdensity in t he IGM, which 
we take to be Aigm = 150 <|Fan et al. iboOrj ). We compute 
/hi (A) at a particular overdensity assuming ionization equi- 
librium, and then integrate over the volume weighted proba- 
bility distribution of the gas density, Pv(A). Since there are 
no measurements of F_i2 at z > 6 we have no independent 
constraints on the mean free path. We therefore adopt an 
upper limit of A m f p = 20 comoving Mpc at z > 6. In prac- 
tice, since Qhii = 1 occurs at z < 6.25 for all three models, 
changing the mean free path (or indeed Pv(A), which is also 
rather uncertain) does not substantially change the (/hi)v 
we derive. As already discussed, the main uncertainties in 
the calculation at z > 6 are the values adopted for a s and 
Chii- 

The evolution of (/hi)v as a function of redshift is 
shown in the left panel of Fig. [9] The constraints on (/hi)v 
at 2 < z < 6, shown by the filled triangles and open di- 
amonds, are computed using the ionization rates inferred 
from the IGM Lya effective optical depth in this work and 
in B05, respectively. When doing this we assume Qhii = 1 
in equation (|25p ; there is strong evidence the h ydrogen in 
the IGM is already hig hly ionized at z < 6 jFan et al. I 
120061 ; iBecker et all 120071 ). The values we infer for the vol- 
ume weighted neutra l fraction a r e in g ood agreement with 
the recent data from iFan et al. I ((2006) at z ~ 6, shown by 
the open circles. The solid line corresponds to (/hi)v com- 



puted at intervals of Az — 0.5 using the updated ionizing 
background model of HM01 for galaxy and quasar emission. 
There is also very good agreement between this model and 
the data based on the IGM Lya effective optical depth. 

Observational constraints on (/hi)v at z > 6 using 
alternative observational techniques are also shown. The 
cross in the left panel of Fig. [9] corresponds to the upper 
limit on (/hi)v at z = 6.3 ob t ained from the spectrum 
of GRB 050904 bv lTotani et al.l l|2006l ). The data point at 
z — 6.25, represented by the inverted t riangle, shows the re- 
cent lo wer limit on (/hi)v measured bv lMesinger fc Haimanl 
|2007l ) using the spectroscopically observed sizes of near- 
zones around quasars combined with modelling of the Gunn- 
Peterson trough damping wing. Note, however, that the es- 
timates of the neutral fraction ba sed on the observed near- 
zone sizes alone ar e ambiguous (iBolton fc Haehneltl 120071 ; 
iMaselli et all 120071 ); the near-zone data is also consistent 
with an IGM which is highly ionized at z ~ 6. The open 
sq uare corresponds t o the upper limit on (/hi)v obtained 
bv lKashikawa et al.l (|2006l ) from the observed evolution in 
the Lycv emitter luminosity function between z = 5.7 and 
2 = 6.5. However, this estimate is also subject to co nsider- 
able uncertainties (|Santos 2004; D iikstra et al.ll2007l ). 

The estimates for (/hi)v computed with equations (|25|l 
and H26|) f° r the three Ni on models at z > 5 are shown 
as the dotted, dashed and dot-dashed lines. The H n filling 
factor, Qhii, is computed assuming Chii = 2. Model 1 is 
inconsistent with all the observational constraints on (/hi)v 
at 2 > 5. Similarly, although model 2 predicts a somewhat 
smaller (/hi)v at z > 5, hydrogen reionization is still not 
complete by z = 6, in disagreement with the neutral fraction 
inferred from the IGM Lya effective optical depth. Hence the 
values of (/hi)v predicted by both model 1 and 2 lie well 
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above the measurements from the Lya data at z = 5 — 6. 
Only model 3 (and also 2b, not shown) is in agreement with 
the Lya data. 

We may also briefly discuss to what extent the three 
models are consistent with current constraints on the inte- 
grated reionization history from the cosmic microwave back- 
ground. The recently re-evaluated value for the electron scat - 
tering optical depth is r c = 0.09 ±0.03 l|Spergel et al.l l2006). 
The electron scattering optical depth may be compute d as 
(e.g. IWvithe fc Loebll2003l : IChoudhurv fc Ferrarall2005h 



= ca T (nH) / dz Qmi(z)(l + z) 3 



(27) 



The cumulative contribution to function of redshift 

is shown in the right hand panel of Fig. [9] for our three 
models. As might be expected, only model 3 is consistent 
with the current constraints on r c . Note again that we have 
arbitrarily assumed that the ionizing emissivity is zero at 
z > 15. Obviously for model 2 and 3 the value of r c would 
be larger for an earlier start of reionization. Free electrons 
left after recombination m ay also increase r e somewhat (see 
IShull fc~Ve nkatcsan 2007 for a recent detailed discussion of 
these issues). 



7.4 A photon-starved extended epoch of 
reionization: Implications for 21cm 
experiments and the photo-evaporation of 
mini-haloes 

Even for a hard spectral index for ionizing photons, a s ~ 1, 
the comoving ionizing emissivity at z = 6 inferred from the 
Lya data is so low that the epoch of reionization has to 
extend over a wide range in redshift. This is obviously ex- 
cellent news for upcoming 21cm experiments as there should 
be plenty of structure observable in 21cm emission over the 
full accessible frequency range of planned experiments. It ap- 
pears likely that a significant fraction of the IGM is still neu- 
tral at 8 = 8, which would bode particularly well for obser- 
vations of H ii regions around early ionizing s ources by up- 
comi n g 21cm experiments such as LOFAR (e .g Wvith e et al.l 
120051 ; IZaroubi fc Silkl 120051 ; iRhook fc Haehnerdl2006h . The 
only downside regarding 21cm experiments is that there is 
little hope for those experiments aimed at a detect ion of a 
global, sharp reionization signal (|Shaver et al.lll999l ). 

An extended photon-starved epoch of reionization also 
has important implications for the photo-evaporation of 
mini-haloes, i.e. dark matter haloes with virial tempera- 
tures below 10 4 K. The virial temperature of these haloes 
is too low for collisional cooling by atomic hydrogen to be 
efficient. Prior to reionization they are thus expected to 
be filled with neutral hydrogen. However, after reionization 
commences they will be photo-evaporated as photoioniza- 
tion raises the gas temperature above the virial temperature 
of the halo. The number of photons per hydrogen atom in 
the halo needed for photo-evaporation, and hence the du- 
ration of the photo-evaporation process itself, will there- 
fore depend on the number of recombinations occurring. 
Illiev et al] l|2005h studied this in detail for a large suite of 
radiative transfer simulations. For haloes at the upper end of 
virial temperatures with total masses of about 1O 7 M0 and 



photo-evaporating ionizing fluxes that corresponds to our 
measured photoionization rate at z = 6, the time needed 
for photo-evaporation approaches a Gyr. Unless the ioniz- 
ing emissivity rises dramatically at z > 6 the photoioniza- 
tion rate towards higher redshift will be even lower. It should 
therefore take until about z = 3 — 4 to fully photo-evaporate 
the most massive mini - haloes . This would support the sug- 
gestion of lAbel fc Mol I 1998) that a significant fraction of 
observed Lyman limit systems, which increase in number 
with increasing redshift, are due to optically thick neutral 
hydrogen in mini-haloes. 



8 SUMMARY AND CONCLUSIONS 

We have obtained new measurements of the metagalactic 
hydrogen ionization rate at z = 5 and 6 using a large suite 
of hydrodynamical simulations combined with rec ent mea- 
surements of the IGM Lya effective optical depth (jSongailal 
l2004l ; lFan et al. 1120061 ). We carefully take into account the 
various systematic errors associated with determining T_i2 
at these redshifts. We find logr H i = -12. 28+0 23 at z = 5 
and logTHi < —12.72 at z = 6, where the largest contri- 
butions to the uncertainties come from the poorly defined 
thermal state of the IGM and measurements of the Lya ef- 
fective optical depth. We also investigate the impact of ion- 
izing background fluctuations on the inferred values of r_i2, 
and find that these add an additional uncertainty of about 
10 per cent. 

Using a physically motivated model for the ionizing pho- 
ton mean free path, we compare our constraints on T_i2 to 
the expected ionization rate calculated from the observed 
galaxy and quasar population at z = 5 and 6. We find that, 
even for conservative estimates regarding the spectral shape 
of the UV emission of these sources, the combined ionizing 
emission from galaxies and quasars is capable of maintain- 
ing the IGM in its highly ionized state if / eS c > 0.2. Ionizing 
emission from star forming galaxies is likely to dominate the 
total ionizing photon budget at these redshifts. The clump- 
ing factor of ionized hydrogen should also be substantially 
less than the often used value of Chii = 30 at z = 5 — 6. 
Our measurements of F_i2 and estimates of A m f p suggest 
Chii < 3 at z — 6, although the clumping factor will increase 
towards lower redshift. 

Using our estimates of the ionizing photon mean free 
path we have turned our measurements of F_i2 from 
the Lya effective optical depth into a measurement of 
the emission rate of ionizing photons per unit comov- 
ing volume at 2 < z < 6. We find that Ni on (z) = 
10 so.5-o.o6( z -6) g -i Mpc -3 ^^1^ Qs = Qb = 3 above 

the Lyman limit frequency. A harder spectral index of 
a s = Qb = 1 would double A^ on . For a s = «b = 3 and 
1 the value of Ni on at z = 6 corresponds to around 1.5 
and 3 photons emitted per hydrogen atom over a time in- 
terval corresponding to the age of the Universe at z = 6, 
respectively. Reionization must therefore have occurred in a 
photon-starved manner unless the ionizing emissivity during 
reionization is substantially larger than at z — 2 — 6. 

We have discussed three simple extrapolations of Ai on 
at z > 6. An ionizing emissivity which decreases rapidly 
towards higher redshift, consistent with rece nt upper limits 
on the UV emissivity from LBGs at z > 6 feouwens et al.1 
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I2005I ), can be ruled out based on the ionization state of the 
IGM at 2 < 6. Consistency with the completion of reion- 
ization before 2 = 6 as inferred from the (/hi)v in the 
IGM, as well as the recently re-determined electron scat- 
tering optical depth, is achievable if the ionizing emissivity 
at 2 > 6 increases above its inferred value at z = 6. Al- 
ternatively, consistency may also be achieved if the 1500A 
luminosity density and the star formation rate remain con- 
stant while / csc increases at 2 > 6, or the spectrum of emit- 
ted ionizing photons becomes increasingly harder at high 
redshift. Recent studies of LBGs at z = 5 suggest that 
typical UV selected star -forming galaxies become younger 
with increasing redshift l|Verma et al.l 12007), which may be 
a hint in this direction. Reionization by an early popula- 
tion of sources such as mini-quasar s or population-Ill stars 
with particularly hard spectra (e.g. IVenkatesan et al.ll200 ' 



ISokasian et alj|2004 iMadau et al.ll2004l ; iRicotti et 81.1120051 ) 
would be an extreme form of such a scenario. 

Assuming e(1500) = 6e L {e.g. MHR99), a con- 
stant ionizing emissivity of around e(1500) ~ 3.8 x 
10 26 (q s /3)(/ csc /0.2)- 1 erg s" 1 Hz -1 Mpc -3 is required for 
reionization before or at z — 6. This is at the upper end of 
some of the recent first attempts to constrain the UV lumi- 
nosit y density at z > 6 feouwens et al.ll2005l; iRichard et al.l 



120061 ; iMannucci et~ai] 120071 ; IStark et al.M2007bl ). There is 
thus reason for optimism regarding future searches for the 
UV emission from high-redshift galaxies. Note, however, 
that the required 1500A luminosity density could be smaller 
if the UV spectrum of the galaxies is harder or the escape 
fraction is larger than we have assumed. 

The photon-starved nature of reionization suggested by 
the Lya forest data at redshift z — 5 and 2 = 6 means 
that the epoch of reionization has to extend over a wide 
redshift range and is unlikely to have been completed much 
before z = 6. In turn this implies a rather early start of 
reionization as suggested independently by the CMB data. 
There should thus be plenty of structure due to the epoch 
of reionization observable in 21cm emission over the full ac- 
cessible frequency range of planned 21cm experiments. The 
photon-starved nature of reionization also means that the 
photo-evaporation of the neutral hydrogen in mini-haloes 
with virial temperatures below 10 4 K should extend to red- 
shifts well below the tail-end of the hydrogen reionization 
epoch. The alternative is that the IGM was rapidly reionized 
at very high redshift by an as yet unidentified population of 
sources which have disappeared by 2 ~ 6. 
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